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1.   Introduction 

A  technique  of  shifting  algorithms  for  tree  partitioning  problems  was  recently  developed 
in  a  sequence  of  papers  ([PS],  [BPS],  [BP],  [ABP])  This  technique  is  a  top-down  greedy 
technique  rather  than  the  bottom  up  approach  used  usually  for  such  problems  (see  e.g., 
[KM] ,  [KH] ). 


An  efficient  implementation  of  the  shifting  algorithm  ([BPS])  for  min-max  tree  partitioning 
is  given.  The  complexity  is  reduced  from  0(Rk  +  kn)  to  0(Rk(k  +  logd)  +  n)  where  a  tree 
of  n  vertices,  radius  of  R  edges,  and  maximum  degree  d  is  partitioned  into  k+ 1   subtrees. 


The  improvement  is  mainly  due  to  a  datj  structure  which  suggests  a  succinct 
representation  for  subsets  of  edges,  of  a  given  tree,  that  preserves  the  interrelation 
between  the  edges  on  the  tree.  Section  2  introduces  both  the  min-max  tree  partitioning 
problem  and  high-level  description  of  the  shifting  algorithm  which  was  given  in  [BPS]. 
Section  3  presents  an  implementation  scheme  of  the  algorithm  which  forms  an 
intermediate- level  description  of  the  algorithm.  Section  4  introduces  two  data-structures. 
It  is  then  shown  how  to  utilize  them  by  the  intermediju  level  description  to  form  a  low- 
level  implementation  of  the  algorithm. 

2.  High-Level  Description 

Let  T(V.E)  be  a  tree  whare  V  is  the  set  of  vertices  and  E  is  the  set  of  (undirected) 
edges.     We  associate  a  nonnegative  weight  w(v)  withi  gvery  vertex  veV.     A  q-partition  of 

the  tree  T  into  q  connected  components  T^,T  T    is  obtained  by  deleting  i,  i<k=q-1  edges 

of  T  where  I5k<n  and  n=|V|.  (The  definition  allows  empty  T,  l^j^q).  The  weight  WfT) 
of  each  component  T^  is  the  sum  of  the  weights  of  its  vertices.  Each  T  is  a  tree,  hence 
we  may  refer  to  it  as  a  subtree  of  T. 

The  min-max  q-partition  problem.     Find  a  q-partition  of  T  minimizing  m^x    WfT) 

Definitions  and  Notations 

For  two  vertices  u.vsv  define  d(u,v),  the  distance  between  u  and  v  in  T,  to  be  the 
number    of    edges    in    the    loop-free    path    from    u    to    v    in    T.  The    radius    R    of    T    is 

min    max       (d(u,v)).  A  vertex  that  yields  the  minimum  in  the  definition  of  R  is  called  a 

{uev>  {v€V> 

center  of  T  Add  an  auxiliary  vertex  r  to  T  with  w(r)  =  0.  Connect  r  to  a  center  by  an 
auxiliary  edge.  Transform  the  new  tree  into  a  rooted  tree  by  choosing  r  to  be  the  root  and 
imposing  a  top-down  direction  on  the  edges.  From  now  on  we  denote  by  T  the  new 
rooted  tree. 

In  this  paper  we  use  the  usual  terminology  of  graph  theory.  If  e  is  a  directed  edge 
incident  from  v^  and  incident  to  v  denoted  by  v  — ^  v  then  we  refer  to  v  as  tail(e) 
and  to  V     as  head(e).     Edge  e  is  said  to  be  the  father  of  edge  e     if  head(e)  =  taiKe  ),  and 


in  this  case,  e.   is  said  to  be  the  son  of  edge  e.  Edges  e.   and  e^  are  said  to  be  brothers 
if  taiKe,)  =  taiKe^.'- 

Edge   e 


path    from    e      to    e_    (resp       e.    to    e.).      Tho   definitions    of    father,    son,    descendent   and 

ancestor  extend  in  the  obvious  way  to  vertices.     Consider  a  function  f:{1.2 k}   — 'E.     Each 

such  function  represents  the  i,  i<k,  edges  L-^ing  deleted  from  T  in  order  to  form  a  (k+1)- 
partition,  where  the  definition  of  a  q-partition  extenas  to  rooted  trees  in  the  obvious  way. 

(The   number    i    is   the   cardinality   of   the   set    {f(1).f(2) f(k)}).      Every   (k+ 1)-partition   of   the 

rooted  tree  T  induces  a  (k+l)-partitiori  on  the'  given  undirected  tree.  Let  f(i)  =  e  for  i, 
1<i<k,  and  edge  e=v^— *"V  ;  then,  we  say  that  cut  i  is  assigned  to  edge  e  and  is  incident 
from  V  We  call  the  function  f  a  configuration  of -nJ!.  (configuration,  in  short).  A 
configuration  g  is  said  to  be  obtainable  from  a  configuration  f  by  a  down-shift  (resp. 
side-shift)  of  a  cut  i,  1<i<k,  if  f(j)  =  g(j)  for  every  j#i.  1<j<k.  and  f(i)  is  a  father  (resp. 
brother)  of  g(i). 

We  further  need  the  notions  of  complete  and  partial  rooted  subtrees:  A  subtree  T  of  T 
is  a  complete  (resp.  partial)  subtref  of  T-.root;ed,at  verte?^  v  Lresp.  edge  v  -*  u)  if  v  is  the 
root  of  T^  and  T^  contains  evsr;  *;on  jlresp.- the^ -son  .^j)  'f  v  together  with  all  the 
descendents  of  these  sons  (resp.  this  sc  .-  -xe       ■   -o     -.CbZ 

-        i      -c   •-)        -   •  •- 

Given  a  configuration  of  cuts  f  let  v— ►  u  .  v— *;  Uj v  j^        be  all  edges  incident  from 

vertex  v  which  are  not  assigned  CL'ts  The  dovvji-comnon.  it  of  v  is  a  subtree  of  the 
complete  subtree  of  T  rooted  at  v;  its  vertices  include  v  and  the  vertices  of  the  down- 
components    of    vertices    u    u  u  .        (Note,    that    this    is    a    recursive    definition).        The 

down-component  of  a  cut  i  is  the  down-component  of  head(f(i)),  and  i  is  called  the  top- 
cut  of  that  component  The  weight  W(i)  of  a  cut  i  is  the  weight  of  its  down  component 
The  down-component  of  an  edge  e  is  the  down-component  of  head(e). 

A  component  T  is  lighter  than  another  component  T  (or  T  is  heavier  than  T)  if  W(J)  < 
W(T). 

The  Shifting  Algorithm 

1.  Initially,  the  k  cuts  are  all  assigned  to  the  edge  incident  from  r. 

2.  While  The  root  component  is  not  a  heaviest  component 
do 

2.1  Down-shift  a  cut  i  of  heaviest  down  component  to  a  vacant  (i.e.  not  assigned 
any  cut)  son-edge  of  a  heaviest  down  component  If  no  such  vacant  edge  exists 
then  halt 

2.2  Let  e  be  the  edge  assigned  to  the  down-shifted  cut  i  before  the  down-shift 


2.3  While  The  eage  e(=v^  — •  v,*  is  vacant  and  \j^*  r 


2.3.1  if  there  is  a  cut  incident  from  v    'such  that  its  down  component  is 
lighter  than  the  down  cornponent  .9f  v, 

then  side-shift  the  lightest-  3uch  cut  tcte   -?:    .-•- 

2.3.2  e*-  the  father  edjd^f  :t:    :".  5l_  :;    2    ^ 

od  '--       ''^    -=-'S    '^-    .'5    =      - 

od 


;2  c£:o-  •    ■^■.•nt-   ens  sJslqmo? 

In  words,  the  algorithm  grhploys  -ai^^)lfefflai^tetoii>  'ifnslructions  2.1,  2.2  and  2.3)  which 
contains  an  internal  loop  (instruction  2.3^'  Gtveri  a  "configuration  of  cuts  let  us  describe  the 
way  it  is  changed  by  one  application  of  the  exteri.-il  loop.  First,  one  of  its  cuts  is  down- 
shifted. Second,  starting  f'-'-  i  the  location  of  this  cut  we  follow  the  path  towards  r  (up 
the  tree)  until  we  first  enct  iter  an  edge  assrjnecT  a  cut  or  arrive  to  r.  Along  this  path 
side-shifts  are  performed.  2&e~:&ik-^cr\ieem  for  "where  to  perform  a  down-shift  or  side- 
shifts  in  the  algorithm  itself.  ?      .  - 

In  [BPS]  it  is  shown  that  the  configuration  which  implies  the  min-max  partition  is 
obtained  before  the  algorithm  halts.     It  is  not  necessarily  the  last  partitioa 

3.  Intermediate -Level  Description 

The    junction-tree    data- structure    which    is    given    in    the    next    section    represents    the 
configurations  obtained  through  the  algorithm.     It  equips  us  with  the  ability  to: 
1.  Update  the  junction  tree  in  0(k)  time  whenever  a  down-shift  occurs 

Z  Update  the  junction  tree  and  compute  in  0(k)  time  all  side-shifts  required  per 
each  down-shift 

Now,  we  describe  additional  information  which  is  kept  and  updated 

For  each  cut  i,  1<i<k,  we  maintain  the  edge  f(i)  (=taii(f(i))  — *  head  (f(i))  which  is  assigned  to 
this  cut  by  the  present  configuration  f  and  the  weight  of  its  down-component  denoted 
W(i). 

We  keep  for  each  vertex  vev,  Z(v)  -  the  weight  of  the  complete  subtree  rooted  at  v, 
and  for  each  edge  esE,  Lie)  -  the  number  of  edges  in  the  directed  path  from  r  to  e  called 


the  level   of   e(1    <  L(e)   <   R^l).      It  takes   0(n)  time  to   compute  Z(v)  and  L(e)   for  the  entire 
tree. 

For  each  veV  we  keep  and  update  a  data-structurs  called  cut-free(v).  It  is  described  in 
the  next  section.  This  data- structure  represents  xhe  edges  e=v— ^■u  such  that  the  partial 
subtree  rooted  at  e  contains  no  cuts,  and  s^fisfius  tl.f    Toll'    ;inc: 

1.  Each  delete  operation  takes  Odog  d  )  time  where  d^  i^  the  number  of  sons  of 
V  in  T. 

2.  Each  insert  operation  takes  0(1)  time. 

3.  It  enables  us  to  get  such  an  edge  e  of  maximum  Z(hccd(e)l  in  0(1)  time. 

It  takes  0(n)  time  to  initialize  this  datn~stfL..;tws'a^:  •■  t;-? . 

In  addition  to  that  we  keep  a  cut-levpl  .arjay.  of  ^ae  k(R+1)  which  is  initially  undefined 
The  entry  (i.h)  will  include  the  h-th  r-iit^f^  jan,  the  dir^^ted  patb,  from  r  to  the  edge  currently 
assigned  to  cut  i.  ^._.^__    ,^  .    ^^..  ^^^ 

Let  us  overview  the  implementation  of  the  sKiffing  algorithnr 


Theorem  V  Each  application  of  the  external  while  \oop  of  the  shifting  algorithm  takes 
O(k+logd)  time,  where  d  is  the  maximum  degree  of  a  vertex  in  T. 

Remark:  We  give  in  this  section  a  proof  which  includes  some  unproved  statements 
related  to  the  junction  tree  and  cut- free  data  structures.  These  statements  are  proved  in 
the  next  section. 

Proof:  There  are  k  cuts.  We  choose  in  0(k)  time  the  cut  i  of  maximum  W(i)  to  be 
down-shifted  in  instruction  2.1.  To  which  edge  should  we  down- shift  cut?  An  edge 
v  — ••v  with  down-component  of  maximum  weight  is  chosen  among  the  candidates 
suggested  by  the  cut-free(v  )  data- structure  (v  =head(f(i))  where  f  is  the  configuration 
before  the  down-shift)  and  candidates  representing  the  partial  subtrees  containing  cuts. 
The  candidate  of  the  first  kind  and  its  weight  are  obtainable  in  0(1)  time.  All  the 
candidates  of  the  second  kind  and  their  weights  are  obtainable  in  0(k)  time.  To  see  this 
we  introduce  an  auxiliary  vector  weight(e)  for  each  eeE  Initialize  weight(e)  to  z(head(e)) 
before  each  applicaiton  of  the  following  procedure.  Upon  termination  of  this  loop  all  the 
entries  of  this  vector  changed  during  the  loop  contain  weights  of  candidates  of  the  second 
kind. 

for  each  j,   1<j<k 
do    et— cut-level(j,Uf(i))+l) 
if  tail(e)=v^ 

then  weight(e)  •>— weight(e)  -  W(j) 
od 

Explanation:     The  cut-level  array  includes  for  each  cut  the  information  whether  it  passed 


through  v,    and  if  yes  through  which  edge  incident  from  v.    it  proceeded 

The  remaining  details  m  order  to  complete  the  proof  that  a  down-shift  takes  Oik)  time  do 
not  require  new  ideas  and  are  therefore  omitted.  Note,  however,  that  if  the  candidate 
suggested  by  the  cut-free(v  )  data-structure  is  chosen  we  need  to  delete  it,  implying 
additional   Odog  d     )   time.      So  ir,Jb-\jdi\(Sa'- 2\t  takes- 0(k   +   log   d    )  time.     The  next  lemma 

deals  with  the  internal  loop  and  thuS  com|iletes  the  proof  of  the  theorem  for  the  external 
loop. 

Lemma  _L  The  total  time  required  by  invocations  of  the  internal  loop  per  each  application 
of  the  external  loop  of  the  aigdri-ltim  is -01k).' 

Proof:  There  are  at  most  0(k)  side-shifts  Tifi'\.^S  internal  loop)  per  each  down-shift  of 
instruction  2.1,  because  each  cut  can  be  side-shifted  at  most  once  as  we  climb  up'  the 
tree  towards  the  root  The  next'sectidh  "jjrefeVWs  ^5n  impiementation  in  which  we  make  up 
to  0(k)  stops  on  our  path'  u^'  the^'trle  to"' cHiBcTc^'for  possible  side-shifts.  In  this 
implementation  the  total  amount  of  time  for  update*  and  computation  of  the  required 
weights  of  down-components  d^es  not  exceed  0(k). .  This  is  due  to  the  second  property 
of  the  junction-tree  data  structure  nign'tiOned  above.  '3oth  the  proof  of  the  theorem  and 
the  lemma  are  completed  at  this  point  prQvjded..,that  we  present  the  promised  properties  of 
the  two  data- structures.  ^    ^^  ....    .         ..-.,     .     , 


4.  Low-level  description 

Given  a  configuration  f  we  define  J(f),  the   junction  tree  of   f  as  follows.     It  is  a  rooted 

tree.       The    set    of    vertices    contains    the    set   of    cuts    {1,2 k}.      In    addition    to    that    it 

contains  vertices  v  e   v  such  that  v  is  a  lowest  common  ancester  of  two  or  more  edges 

of  the   set    {f(1),   f(2) f(k)}.     The  latter  are  called   junction  vertices  while  the  former  are 

called  cut  vertices.  There  is  an  edge  from  a  cut  vertex  i  to  a  cut  vertex  j  (resp.  junction 
vertex  v)  if  there  is  a  directed  path  from  f(i)  to  f(j)  (resp.  v)  in  T  that  does  not  include  any 
other  vertex  of  J(f).  There  is  an  edge  from  a  junction  vertex  v  to  a  cut  vertex  i  (resp.  a 
junction  vertex  u)  if  there  is  a  directed  path  from  v  to  f(i)  (resp.  u)  in  T  that  does  not 
include  any  other  -vertex  of  J(f).  Figure  1(a)  describes  an  example  of  T  and  f.  Figure  Kb) 
gives  J(f).  Note  that  a  cut  vertex  has  at  most  one  son  and  a  junction  tree  has  at  most 
2k- 1   vertices. 

Let  us  show  that  any  down-shift  or  all  side-shifts  per  down  shift  can  be  reflected  on 
the  junction  tree  data  structure  in  0(k)  time. 

Below,  we  use  the  convention  that  f(resp.  g)  is  the  configuration  of  cuts  before  (resp. 
after)  the  described  shift  The  letters  v  (resp.  i.  resp.  a.  ^)  represent  junction  vertices 
(resp.  cut  vertices  representing  cuts  being  shifted,  resp.  portions  of  a  junction  tree)  of  J(f) 
or  J(g).  Let  us  describe  possible  changes  in  J(f)  following  one  down-shift  Figure  2 
describes  the  case  where  v=head(f(i))  is  a  junction  vertex  (see  Figure  2(a)).  There  are  two 
possibilities:      either   i   is   down-shifted  to   an   edge  of   T  on  the   path   from  v  to   some  cut 


(Figure  2(b))  or  not  (Figure  2(c)).  Figure  3  describes  the  case  where  head  (f(i))  is  not  a 
junction  vertex  (see  Figure  3(a))  Figure  3(b)  describes,  actually  two  possibilities:  the 
junction  tree  portion  H  is  empty  or  it  is  not  empty  and  cut  i  was  down-shifted  to  the  next 
edge  in  T  on  the  path  from  f(i)  to  the  (vertex  or  edge)  in  T  represented  by  the  son  of  i  in 
J(f).  Figure  3(c)  describes  the  other  case  of  Figure  3.  Note  that  in  the  last  case 
v=head(f(i))  becomes  a  junction  vertex. 

How  to  actually  perform  each  such  change  in  (jik)  time?  Notice  that  in  the  cases 
described  in  Figures  2(b)  and  3(b)  (when  ^  is  not  e.ipty)  the  edge  to  which  the  cut  is 
down-shifted  is  not  taken  from  the  cut-free  data-s+'MCture.  The  computation  of  weights 
of  the  candidates  required  for  these  cases  can  be  slightly  modified  to  yield  the  edge  in  J(f) 
to  which  cut-vertex  i  is  moved:  start  from  a  cut  vertex  that  changed  weight  (g(i)-  the 
down  shifted  cut)  in  the  implementation  of  instruction  2. 1  described  in  the  proof  of 
Theorem  1.  Follow  the  path  up  J(f)  in  order 'W  identify  the  edge  of  J(f)  to  which  cut 
vertex  i  is  moved  to  form  J(g).  The  otliP''  changes  required  in  this  case  are  simple.  In  all 
the  other  cases  of  down-shifts  cut  vi..,..^  i  'creates'  a  new  edge  in  J(g)  which  is  trivial  to 
implement 

The  possible  changes  in  J(f)  foll'"ving  a  sido-shift  are  described  in  figures  4  and  5. 
Figure  4  describes  the  case  where  cut  vortex  i  does  not  have  descendents  in  J(f)  and 
junction  vertex  v=tail(f(i))  has  only  two  sons  (including  i).  This  causes  the  elimination  of 
junction  vertex  v  in  J(g),  see  Figure  4(b).  Fipr-e  5  includes  the  other  cases:  at  least  one 
of  the  junction  tree  portions  a_  and  a-  is  :.«.t  empty.  The  result  is  seen  in  Figure  5(b). 
The  implementation  here  is  trivial  assuming  thi.t  we  know  which  side-shift  to  perform.  We 
see  below  that  we  use  the  junction-tree  structure  to  show  us  which  side-shifts  to 
perform. 

Implementation  of  the  internal  loop 

We  keep  an  auxiliary  variable  W  which  is  set  to  zero  before  each  application  of 
instruction  2.2.  Instead  of  climbing  on  T  we  climb  on  the  junction  tree.  At  each  point, 
before  climbing  further  on  the  junction  tree,  W  is  set  to  equal  the  sum  of  weights  of  all 
down  components  of  cuts  which  are  assigned  to  edge-descendents  of  our  present 
position  in  the  climbing.  We  initialize  W  at  the  down-shifted  cut  For  instance,  in  the  case 
described  in  Figure  2(a)  W  is  initially  the  sum  of  alt  weights  of  down-components  of  cuts 

that    belong    to    the    junction    tree    portions    a^,    a  a.    and    the    weight    of    the    down 

component  of  cut  i.  The  implementation  of  the  internal  loop  goes  then  as  follows. 

Let  vertex  u  denote  our  present  position  in  the  climbing.  Vertex  u  belongs  to  the 
junction-tree  data- structure.  Initially  u  is  the  father,  in  the  junction-tree,  of  the  last  down- 
shifted cut 

While  u  is  a  junction  vertex 

do 
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(let  u — 't  be  the  edge  in  the  junction  tree  on  which  we  climbed  to  u,  let  u^v  be 
the  first  edge  in  the  path  in  T,  which  corresponds  to  u— t  in  the  junction  tree.  The 
weight  of  the  down-component  rooted  at  v  is  Z(u)-W) 

if  there  is  a  cut  vertex  j  which  is  a  son  of  u  in  the  junction  tree.  f(j)  is  incident 
from  u  in  T  and  the  weight  of  cut  j  is  lighter  than  Z(v)-W 

then  the  lightest  such  cut  is  chosen  and  it  is  side-shifted  to  u — ►v,  and  W-«— Z(u)  + 
J]W(j).  (The  sum  is  over  cuts  j  which  were  descendents  of  u,  but  not  v,  before  the 
last  side-shift,  not  including  th^j  cut  which  was  side-shifted  above) 

else  W-«— W+EW(j)  (The  sum  is  over  cuts  j  which  are  descendent  of  u,  but  not  v). 

If  u  is  the  root  of  the  junction  tree 

then  go  to  the  next  performance  of  the  external  loop 

else  u^—  the  father  of  u  in  the  junction  tree 
od 

(we  arrive  at  a  cut  vertex  i)  W(iH— Z(head  (f(i)))  -  W. 
Some  remarks  have  to  be  added 

1.  The  connection  between  locations  in  the  junction  tree  and  their  corresponding 
locations  in  T  is  done  by  using  at  any  time  the  level  L  in  conjunction  with  the  row  of  the 
cut-level  array  corresponding  to  the  cut  which  was  down-shifted 

2.  Since  we  never  go  back  to  a  portion  of  the  junction-tree  which  was  vistted  before, 
the  computation  time  is  bounded  by  time  which  is  proportional  to  the  size  of  the  junction 
tree-O(k).  Note  that  during  the  computation  of  W  we  only  visit  new  edges  of  the  junction 
tree  (due  to  the  reservation;   descendents  of  u,  but  not  v'). 

The  cut-free  data- structure 

The  cut-free(v)  (for  every  v  e  V)  data- structure  is  a  priority  queue.  It  is  actually  a 
combination  of  two  standard  data-structures:  a  heap  and  a  stack.  Initially  the  heap  part  of 
cut-free(v)  contains  all  the  edges  v— *u  in  T.  The  heap  is  organized  according  to  Z(u)  with 
the  son  edge  of  heaviest  Z{u)  at  the  top.  It  takes  0(n)  operations  to  initialize  the  heaps,  of 
all  vertices  in  T. 

To  remind  the  reader,  the  cut-free(v)  data-structure  represents  the  edges  v-*u  in  T  for 
which   the   partial   subtree  rooted   at  v— 'u   contains   no   cuts.      Our   heap   contains   initially   d 


elements  where  d     is  the  number  of   sons  of   v,   VT^r   in  T    The  deletion  of  an  element  from 
the  heap  and  the  update  that  follows  takes  Odg  d  )  time. 

We  still  have  to  consider  the  possibility  of  returning  an  edge  v— ►u  to  cut-free(v)  while  a 
side-shift  of  a  cut  from  (v,u)  is  performed  and  the  partial  subtree  rooted  at  v— *u  contains 
no  cuts.  We  add  a  stack  to  the  cut-free  data-structure  to  contain  such  edges.  Say  that  a 
down-shift  of  a  cut  i  on  the  edge  w— ♦•v  is  decided,  and  a  choice  to  which  edge  to  down- 
shift should  be  made.  For  this  we  have  to  consider  the  candidate  of  the  stack  as  well  as 
the  candidate  of  the  heap  of  cut-free(v).  Any  access  to  a  stack  for  insertion  or  deletion 
purposes  require  0(1).  The  following  lemma  shows  that  the  top  of  the  stack  of  cut-free(v) 
contains  a  heaviest  partial  subtree. 

Lemma  2:  The  edges  v— nj  in  the  stack  of  cut-free(v)  are  stored  in  non-decreasing  order 
of  the  weights  Z(u)  of  the  partial  subtrees  which  are  rooted  at  them. 

Proof.  Suppose  a  cut  i  was  side-shifted  from  an  edge  v— *u  which  was  then  inserted 
into  the  stack  and  later  another  cut  i  was  side-shifted  from  an  edge  v— mj  which  was 
then  inserted  into  the  stack.  If  i  was  at  the  edge  v— *u-  when  i^  was  side-shifted  from 
V— nj  then  by  the  shifting  algorithm  2(u  )  £  2(u_).  Otherwise,  i_  was  down-shifted  to 
V— ♦u.  while  V— *u  was  vacant  Again  by  the  shifting  algorithm  2(u  )  <  Z(u_).  Hence  the 
order  of  storing  edges  in  the  stack  satisfies  the  Lemma     |    | 

The  number  of  down-shift  operations  of  the  external  loop  is  bounded  by  kR.  Therefore 
Theorem  1  implies  that  the  running  time  of  the  present  implementation  of  the  min-max 
shifting  algorithm  is  O(n+Rk(k+Iog  d)).    ' 
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